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Abstract. The Sznajd model, which describes opinion formation and social influence, is treated analytically 
on a complete graph. We prove the existence of the phase transition in the original formulation of the 
model, while for the Ochrombel modification we find smooth behaviour without transition. We calculate 
the average time to reach the stationary state as well as the exponential tail of its probability distribution. 
An analytical argument for the observed 1/n dependence in the distribution of votes in Brazilian elections 
is provided. 

PACS. 89.65.-s Social and economic systems - 05. 40. -a Fluctuation phenomena, random processes, noise, 
and Brownian motion - 02.50.-r Probability theory, stochastic processes, and statistics 



1 Introduction 

There is significant convergence between statistical physics 
and mathematical sociology in approaches to their respec- 
tive fields p[]. Ising model, the single most studied statis- 
tical physics model, finds its numerous applications in so- 
ciophysics simulations. Conversely, sociologically inspired 
models pose new challenges to statistical physics. We be- 
lieve this is the case of the Sznajd model we arc studying 
here. 

The model of K. Sznajd- Weron and J. Sznajd [2] was 
designed to explain certain features of opinion dynamics. 
The slogan "United we stand, divided we fall" lead to 
simple dynamics, in which individuals placed on a lattice 
(one-dimensional in the first version) can choose between 
two opinions (political parties, products etc.) and in each 
update step a pair of neighbours sharing common opinion 
persuade their neighbours to join their opinion. Therefore, 
it was noted that contrary to the Ising or voter [3] models, 
information does not flow from the neighbourhood to the 
selected spin, but conversely, it flows out from the selected 
cluster to its neighbours. 

The model initiated a surge of immediate interest 

and the results of numerical simulations can be briefly 
summarised as follows. The results do not depend much 
on the spatial dimensionality or on the type of the neigh- 
bourhood selected . In the case of q choices of opinion, 
the system has q obvious homogeneous stationary (absorb- 
ing) states, where all individuals choose the same opinion. 
There is no way to go out of the homogeneous state, so it is 



an attractor of the dynamics. This is reminiscent of a zero- 
temperature dynamics, which in Ising model leads to rich 
behaviour |2fi| . However, in the Sznajd model, the pos- 
sible metastable states, like the "antiferromagnetic" con- 
figuration have negligible probability to occur, unless we 
introduce explicitly also an "antiferromagnetic" dynamic 
rule as it was used in the very first formulation [S]. 

The case q — 2 was studied mostly, denoting the opin- 
ions by Ising variables +1 and —1. The probability of hit- 
ting the stationary state of all +1 (or, complementary, 
all —1) was studied, depending on the initial fraction p 
of the individuals choosing +1. Sharp transition was ob- 
served at value p = 0.5 JTj; for p > 0.5 the probability 
to reach eventually the state of all opinions +1 is close to 
one, while for p < 0.5 it is negligible, which can be inter- 
preted as a dynamical phase transition. The distribution 
of times needed to reach the stationary state was mea- 
sured, revealing a peak followed by relatively fast decay. 
This means that the average hitting time is a well-defined 
quantity 

It was also found in one and two-dimensional lattices 
that the fraction of individuals who never changed opin- 
ion decays as a power with time, similarly to Ising model. 
While the exponent in one dimension agrees with the Ising 
case, the two-dimensional Sznajd model gives different ex- 
ponent than Ising model, indicating different dynamical 
universality class 1 1 3| . Also the waiting time between two 
subsequent opinion changes is distributed according to a 
power- law P]. 



e-mail: slanina@fzu.cz 



Among other studies, let us mention the influence of 
advertising effects |18U19| and price formation [20] . Long- 
range interactions were studied in |21j . 
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In a very short but intriguing note |22| Ochrombcl sug- 
gested a drastic simplification of the Sznajd model. In the 
Ochrombel version it is not necessary to have a cluster of 
identical opinions. Any individual is capable to convince 
her neighbours to select the same opinion. This model was 
reported to share all essential features of the original Sz- 
najd model, only the phase transition in the probability 
of hitting the state of all +1 at p = 0.5 is absent. 

The Sznajd model was also used to model the election 
process. There is recent empirical evidence from Brazilian 
elections |27U28ll2T?] that the distribution of votes per can- 
didate follows a power-law, more specifically P(n) ~ 1/n, 
where n is the number of votes. This result was reproduced 
in a study rj] based on Sznajd model on a scale-free net- 
work HEM- 

The dynamics of elections was thoroughly investigated 
by S. Galam 33 34 35 36 , showing that majority rule ap- 
plied on sufficiently many hierarchical levels leads to a 
homogeneous "totalitarian" state with one opinion per- 
vading the whole system. 

Other approaches to physical modelling of opinion dy- 
namics were also investigated |37U38j and among them es- 
pecially the Axelrod model, which was found to have rich 
behaviour from the statistical physics point of view |39l 

EEll 

We should also mention the well studied voter model 
3,42,43,44 , which is very similar in spirit to the Sznajd 
model. Indeed, the relation of the two models was studied 
e. g. in and it seems that Sznajd model reduces to the 
voter model at least for certain setups (especially using 
the Ochrombcl simplification on a complete graph) while 
for others the voter model can be generalised so that it 
includes the rules of Sznajd model as a special case. In fact, 
similar analysis to that presented here was performed for 
voter model, contact process and related processes in |4fi) . 
The persistence properties of the voter model on complete 
graph were studied in 02] ■ 

Very recently a "Majority rule" model, sharing some 
features with Sznajd model, was introduced and studied in 
|47] and its generalisation to the Majority-Minority model 
|48j gives in the mean-field approximation results closely 
related to ours. 



2 Formulation of the model and its 
simplifications 

2.1 General scheme 

In the original formulation of the Sznajd model, the "united 
we stand" principle is often stressed [2"llll| . It means that 
only a cluster of identical opinions can spread the same 
opinion toward its neighbours. However, this principle was 
relaxed in the Ochrombcl simplification [22| without qual- 
itatively affecting many of the results (except the presence 
of the phase transition) . We will propose some other sim- 
plifications here, supposing the results remain robust. 

Let us have N agents, each of which can be in one of 
q states (opinions) a G S. We may for example think of a 



g-state Potts model variables. Each agent sits on a node 
of a social network, and they can interact along the edges 
with their nearest neighbours. 

The opinion of the agent i is denoted (Tj. The state of 
the system is described by the set of opinions of all the 
agents, S = [<r 1 ,a 2 , ~.,<tn]- 

The variable £(t) performs a discrete-time Markov 
process, whose transition probabilities from time t to t + 1 
differ in various cases, which will be specified in the fol- 
lowing. 



2.2 Case I: two against one 

The first case investigated, which we will sometimes call 
"two against one" , generalises and simultaneously simpli- 
fies the various versions introduced in [TJ . The main dif- 
ference is in the fact that we will change at maximum one 
agent at each time step. This may not significantly change 
the behaviour, as the various choices of neighbourhood in 
[TTj exhibit only little difference. 

Our algorithm will iterate the following three steps. 
First, choose randomly an agent i. Then, choose randomly 
one of its neighbours, say j. If <Ji(t) ^ o~j(t), nothing hap- 
pens. However, if o~i(t) = o~j(t), we will choose randomly 
one of the common neighbours of both i and j, say fc, and 
set ak(t + 1) = (Ji(t). We may also write it schematically 
as reactions AAB -> AAA, BBA -> BBB. 



2.3 Case II: Ochrombel simplification 

In this case, we do not need to have two neighbours in 
the same state. Everybody can influence each of its neigh- 
bours. We choose an agent i at random. Then, choose j 
randomly among neighbours and set o~j(t + 1) = Ui(t). 
Therefore, the process may be written as AB — > AA, 
BA — > BB. In fact, on fully connected network the Ochro- 
mbcl simplification is equivalent to voter model, whose 
dynamical properties were studied e. g. in 12 . 

As an obvious observation we can note that both in 
case I and case II the uniform states, with all (Tj equal, 
are stable under the dynamics. However, we can expect 
variety of metastable states in the case I, in which there 
are no pairs of neighbours in the same state, therefore the 
dynamics does not proceed any further. 



3 On a fully-connected network 

We will approximate the complex social network by the 
fully-connected network (the complete graph) of N nodes. 
Here, any two agents are neighbours; in the case I we sim- 
ply choose three agents k at random and in the case 
II two agents i,j at random. Note that the order in which 
they arc chosen matters. This makes our process differ- 
ent e. g. from the majority |47| or majority-minority |48| 
models, although on fully connected network the differ- 
ence may consist only in rescaling certain variables. 
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We will call this setup a mean-field approximation in 
the same sense as the Ising model on the complete graph 
can be considered as an approximation for Ising model on 
hypercubic lattice of high dimensionality. Of course this is 
not a good approximation to the original one-dimensional 
formulation of the Sznajd model [2], but we believe it 
is appropriate for much more realistic studies of Sznajd 
model on complex networks [HlI6ll23j . We refer the reader 
to | Appendix A| for a more formal definition of the Sznajd 
model on an arbitrary graph. 

In fully-connected network the state of the system is 
fully described by the occupation numbers N a = Yli=i <5a>o-, 
or cquivalcntly the densities n„ = N a /N, for each opin- 
ion a G S. The dynamics of these occupation numbers 
fully describes the evolution of the system. As the total 
number of nodes is conserved, there are q — 1 independent 
dynamical variables. 

Let us start with the case II (Ochrombel simplification) 
with only two opinions, q = 2. The variable a can assume 
only two values, denoted a = ±1 for convenience. Indeed, 
we are effectively working with Ising spins. The state is 
described by one dynamical variable only, which will be 
taken as a "magnetisation" , 



to 



N+-N- 
N 



(1) 



In one step of the dynamics, three events can happen. 
The magnetisation may remain constant or it can change 
by ±2/N. The probabilities of these three events can be 
easily calculated 



(2) 



Prob | m ^m+|| = i(l-m 2 ) (l + j^j 
Prob {to^to-|} = 1(I-to 2 ) (1 + ^ 
Prob{m^m}= Q (l + ™ 2 ) - 1) (1 + -L. 



Our objective is writing the master equation for the 
probability density of the random variable m(t), which we 
denote P m . It can be found easily in the thermodynamic 
limit N — > 00. Indeed, we find that the time should be 
rescaled as 

t = N 2 r (3) 

in the thermodynamic limit. Then the probability density 
evolves according to the partial differential equation 



d_ 



P m (m,T) 



di 



[(l-m 2 )P m (m,r)] 



(4) 



The latter equation describes in principle fully the evolu- 
tion of the Sznajd model in Ochrombel simplification on 
a complete graph. It has the form of a diffusion equation 
with position-dependent diffusion constant. 

Let us turn now to the case I (original Sznajd model), 
again with q = 2. Wc may repeat step by step the consider- 
ations made above for the case II. Namely, our dynamical 
variable will be again the magnetisation to which may ei- 
ther remain unchanged or change by ±2/N in one step. 



For the probabilities of these events we can find formulae 
analogous to @ 



2 1 1-to 2 / I+3m 

2 1 (I - m 2 ) / 

Prob < m — * m > = 1 — to + 

Nj 8 V 

I — 3 m 



N 



(5) 



(1 - TO 2 ) / 1 

Prob{TO -h. m} = 1 - V - ' 1 + - 



N 

where the terms of order 1/N 2 are neglected. Note that 
the probabilities of changes ±2/N are not symmetric, con- 
trary to the previous case (II). This fact has all-important 
consequences. We will see later that it is responsible for the 
fact that the original Sznajd model exhibits phase tran- 
sition, while in Ochrombel simplification the transition is 
absent. 

A more immediate consequence is that the time must 
be rescaled differently, in order to get sensible thermody- 
namic limit, namely 



t = 2N t 



(6) 



The second consequence is that the equation for P m (m, r) 
contains first derivative with respect to to, representing a 
pure drift in magnetisation: 

— P m (m,T) = -—[(1- m 2 )mP m (m,T)] . (7) 

Contrary to the previous case Q the diffusion term, con- 
taining the second derivative in m, represents only the 
finite-size correction to the drift term. However, this cor- 
rection may dominate close to points to = ±1 and to = 
where the drift velocity becomes zero. 

Next case investigated will be the case II with arbitrary 
value of q. Moreover, we will assume that the number of 
opinions is large, q 1. Let us define the distribution of 
occupation numbers 



N 

D(n) = — y 5 (n — n a ) 

9 £ri 



(8) 



where 5(x) = I for x = and zero elsewhere. It would be 
much more difficult to write the full dynamic equation for 
D(n). Therefore, we use the approximation which replaces 
the distribution D(n) by its configuration average P n (n) = 
(D(n)). In the limit N — > 00 and q — > 00 and substituting 
the variable x = 2n — 1 we arrive at the equation 

JL Pn ( X ,T) = -^[(l-X 2 )P n (x,T)] . (9) 

The time is rescaled again according to the Eq. JSJ). We 
can see that the equations Q and have identical form, 
although the interpretation of variables is different. We 
can therefore solve the two cases simultaneously. This will 
be performed in the next section. 
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4 Solution of the dynamics 
4.1 Two against one: case I 

The case I, q = 2 is described by the equation 

9 P(x,t) = --?-[(1-x 2 )xP(x,t)] . (10) 



ch 



It can be easily verified that the solution has the following 
general form 



P(x,r) = [(l-x 2 )x]- 1 f(e 



r) (11) 



for arbitrary function f(y). The form of the function f(y) 
is given by initial conditions. For example if the initial 
condition is a (5-function, it keeps the same form during 
the evolution, only the location shifts in time. This way 
we could in principle calculate, how long it takes to reach 
the edges of the interval from given initial position. This 
would be the time to reach the stationary state. However, 
it comes out that the time needed blows up. The reason 
comes from the infinite-size limit N — > oo. Indeed, very 
close to the points x = ±1 the finite-size effects take over. 

We can estimate the average time needed to reach the 
stationary state in finite system by the following consider- 
ation. In fact, the equation l|10[) describes the drift which 
pushes the system toward the stationary state, but ne- 
glects the effect of diffusion, which becomes important 
at a distance ~ 1/N from the points x = ±1. There- 
fore, we must calculate the time necessary for the drift 
to drive the system to the point ±(1 — 1/N). The initial 
fraction p of opinions +1 corresponds to the initial condi- 
tion xq = 2p — 1 and from the formula jTTJ) we have the 
following estimate for the average time (r s t) to reach the 
stationary state 

It is also possible to include the correction terms of 
order 0(1/N) into Eq. 1|10|> and deduce the equation for 
the average time to reach the absorbing state (T s t)(xo) on 
condition that the process started at initial position xq. 
Following the general scheme we obtain a second-order 
ordinary differential equation 



3 \ d 
+ I(l-,o 2 )^(r st )(,o) = -l. 



The solution of (|T3J) is 

(r s t) (xo) = N 



(13) 



x r O e ^p-z 2 



-1 Jy 



1-z 2 



dze~^ y 'dy. (14) 



Indeed, for xq not too close to either of the points xo = 
— 1, 0, 1 (the distance must be large compared to 1/N) we 



obtain from the formula 1)14(1 an approximate expression of 
the form given in l|12[> . Another way to obtain the same p 
dependence as in (|12(l is to omit the 0(1/N) terms in the 
equation fT3*jl and solve the first-order differential equa- 
tion. In this case, however, we lose any information about 
the dependence on N. We should also note that a result 
essentially equivalent to Eq. ((12JI was obtained also in 02j ■ 

It is rather interesting to observe that the determinis- 
tic dynamics of Galam model 34,36 leads to a formula 
very similar to (|12[) . while the interpretation of the time 
variable is totally different: in Galam model it represents 
the number of hierarchical levels on which the majority 
rule is iterated. 

It would be desirable to calculate the full probability 
distribution for the time to reach the stationary state r st 
and not only the average. That is possible using again 
the formalism of adjoint equation when we introduce 
the 1/N corrections to Eq. i|10|) but the resulting partial 
differential equation is difficult to solve explicitly. Instead, 
we estimate the exponential tail of the distribution by a 
simple consideration. 

Indeed, after the drift had pushed the system to the 
state in which there is only single spin —1 immersed in 
a sea of all +l-s it finally comes into uniform stationary 
state if the first pair of spins chosen is both +1 and the 
third one is the single —1. This choice has probability 
~ 1/N. Therefore, the relaxation time toward the uniform 
state is i re i a x — N and using the scaling JBJ) we have for 
the tail of the distribution 



P(r at ) ~ exp(- 



Tst 
Trelax 



), Tst 



with 



Trcla 



(15) 



(16) 



The most important observation we can draw from the 
solution i|ll|) is the presence of the dynamic phase transi- 
tion, as observed in numerical simulations. Indeed, start- 
ing with any fixed positive magnetisation, we have initial 
condition P(x,0) = 5(x — Xo), xq > 0, and the drift ex- 
pressed by Eq. (|11|) always take us to the state with all 
agents having opinion +1, while from any state with neg- 
ative magnetisation the drift leads the system eventually 
to the state with all agents having opinion —1 and the 
probability of ending in the state of all +1 is therefore 
P+ = 9(p — 1/2). The possible deviations from this rule 
close to the zero magnetisation (i. c. p = 0.5) are due 
to the finite size effects, which are neglected in l|10fl . The 
presence of the phase transition is also indicated by the di- 
vergence of the average time to reach the stationary state 
(O for p -► 1/2. 



4.2 Ochrombel simplification: case II 

The equation 

§ ? P M = ^[(l-x 2 )P(x,r)] (17) 
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describes both the case II, g = 2 and II, q ^> 1, only the 
interpretation of the variable x differ: in the former case it 
corresponds to the magnetisation, while in the latter case 
it is shifted percentage of votes. By solving Eq. (|17|1 we 
treat simultaneously both cases. 

The equation of the form 117|) was already studied in 
variety of contexts, e. g. population genetics |5()U51| or re- 
action kinetics |52| and can be tackled by standard meth- 
ods developed for Fokker-Planck equation. 

Indeed, we look for the solution using the expansion in 
eigenvectors. We can write (II 71) it in the form -^pP(x, r) = 
£P(x,t) where the linear operator C acts as (£/)(x) = 

[(1 — x 2 ) f(x)] . We therefore need to find the set of 
eigenvectors of C. Denoting <P c (x) the eigenvector corre- 
sponding to the eigenvalue — c, we have the following equa- 
tion 

{l-x 2 )^{x)-Ax$' c {x) + {c-2)<P c {x) = 0. (18) 
The full solution of (|17|l can be then expanded as 

P(i,r) = ^4 c e- CT # c (x) (19) 

c 

with coefficients A c determined from the initial condition. 

Important question to be settled prior to the attempt 
for solution is, what is the appropriate space of functions 
<P(x). First, the interpretation of these functions as prob- 
ability densities sets the requirement that it must be nor- 
malisablc: J <P(x) dx < oo. Second, only the interval x G 
[—1,1] is relevant, so <P(x) = outside this interval. Fi- 
nally, we should anticipate the possibility that (5-functions 
appear in the solution, namely located at x = ±1, because 
the uniform states, with all sites carrying the same spin 
value, are stable under the dynamics. 

We therefore look for the solution of (|18|) in the space of 
distributions (i.e. linear functionals on sufficiently differ- 
entiable functions) with support restricted to the interval 
[-1,1]- 

It is straightforward to find the eigenvectors corre- 
sponding to eigenvalue c = 0, i. e. the stationary so- 
lutions of the equation (|17|l . They arc composed of 8- 
functions only. In fact, the corresponding cigensubspace 
is two-dimensional and the base vectors can be chosen as 

$oi=6(x-l) , $ Q2 = 5(x + l) . (20) 

For c/Owe first decompose the solution in ordinary 
function of x plus a pair of 5-functions, namely 

<P c = <p c+ 6(x-l) + (f> c -5(x + l)+ 

+ 4> c (x)0(x-l)0(x+l) ( ' 

where C+ and C _ are real numbers and 4> c {x) is a real 
doubly diffcrentiablc function. Then, the equation l|18|) 
translates into equation for (f> c (x) 

(1 - x 2 ) 4>"{x) - Ax <j)' c {x) + (c - 2) 4> c (x) = (22) 

accompanied by two other conditions 

c 

Hm 4>c{x) = --z<f> c ± • (23) 

x— >±1 Z 



The general solution of the equation (|22J) exhibits be- 
haviour 4> c (x) ~ (l=Fa;) a ata;^±l, where either a = or 
a = —1. However, the latter case should be excluded, as 
it gives non-normalisable probability distribution. In fact 
it is the condition of normalisability that determines all 
possible eigenvalues c. The solution of (|22|l with correct 
behaviour at x — > ±1 can be expressed in Gegenbauer 
polynomials 52,53,54 . The eigenvalues are c = cj = 
(I + + 2) for I = 0, 1,2, ... . An elementary solution 
and the table of several lowest polynomials is presented in 
|Appcndix B| 

It is important to note that for any eigenvalue c > 
we have 

J<P c {x)dx = Jx<P c (x)dx = . (24) 

The consequence is that both J P(x,r)dx and J xP(x,r)dx 
are independent of time. While the first conservation law 
expresses simply the conservation of probability, the sec- 
ond one is a non-trivial consequence of the model dy- 
namics. Mathematically it is related to the fact that the 
eigenspace corresponding to zero eigenvalue is two-dime- 
nsional. 

Thus, we found the set of right eigenvectors of the op- 
erator £. For practical solution we still need to estab- 
lish the coefficients A c in Eq. I|19(l . To this end we need 
also the set of left eigenvectors of C, checking simulta- 
neously that the set of left and right eigenvalues coin- 
cide. First, we need to establish the adjoint operator to 
£, defined by usual relation (£f\g) = (f\C T g). While C 
acts on the space of distributions, its adjoint C T acts on 
the corresponding dual space, which is the space of suf- 
ficiently diffcrentiablc functions. Straightforward algebra 
gives (£ T g)(x) = (1 — x 2 ) g"{x) which implies the follow- 
ing equation for the left eigenvectors 

(1 -x 2 ) il>'l(x) +c^ c (x) = . (25) 

We find again that for c = the eigensubspace is two- 
dimensional. We can choose the basis vectors so that they 
are mutually ortho-normal to the pair of right eigenvectors 
(l2T)jl . namely 

1&Oi = ^(l + aO , ^02 = ^(1- x). (26) 

The solutions of l|25|) for c > with proper boundary 
conditions are again polynomials presented in more detail 
in |Appcndix B| 

The coefficients in the solution ()19(l with initial condi- 
tion P(x,0) = are then calculated as 

A = J Pq(x)Mx) dx 
J 4>c(x)4> c (x) dx 

From the solution (|19|) we can deduce an important 
feature for the distribution of waiting times needed to 
reach the stationary state. Indeed, if P s t(i~) is the proba- 
bility density for ending at time r in the stationary frozen 
configuration with all agents in the same state, we can 
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express the probability that the stationary configuration 
was not reached before time r as 



lim 



-1+e 



(28) 



P(x, r)dx 



We can see that only the <5-function components of the 
eigenvectors <P c (x) in the expansion l|19(l contribute to 
P s ^(r st ). More explicitly we find 



P>(t)=^2A c 



,(-!) + &(i) c _ c 



(29) 



A — 



0.01 



0.001 



c>0 




As the spectrum of eigenvalues is discrete, for long times 
only the lowest non-zero c (equal to cq = 2) is relevant. 
Therefore, the distribution of waiting times will have an 
exponential tail P^{t) ~ c~ 2r , r — > oo. For initial con- 
dition Po(x) = 5(x — xq) we can easily compute also the 
prefactor in the leading term for large r. Indeed, from lt2"7|l 
we get Ai and finally obtain 



Plkr)^%{\ 



xl)e- 2 \ r 



oo 



(30) 



As the functions 4> c {x) are odd for c = c\ with odd I, we 
should expect that the corrections to the formula (|30|l will 
be governed by the second next eigenvalue c-i = 12. We will 
see later how it can be checked in numerical simulations. 

As in the case I the average time (r st ) (xq) to reach the 
absorbing state when starting at position Xq can be ob- 
tained, using the general formalism |49j . from the equation 



which can be solved easily 



/ \, \ x l + x 11 

2 1 — Xq 2 



XQ 



(31) 



(32) 



Fig. 1. Probability of reaching the stationary state in time 
larger than r, for case I, q = 2, iV = 2000. The values of initial 
fraction p of opinions +1 are 0.1 (+) 0.2 (x) and 0.7 (0). 



(see also jHSES])- The method of adjoint equation p!??ll53"| 
can be used to calculate the distribution of times to reach 
the absorbing state, when starting from initial position 
at x = Xq, yielding results equivalent to our direct cal- 
culation. Indeed, inserting the initial condition Po(x) = 
S(x — xo) into Q27JI we can see that the expression l|29ll 
represents an expansion in the eigenvectors ip c (xo) of the 
adjoint operator C T taken at point xq- 

Contrary to the case I, we do not observe any phase 
transition here. This is due to the conservation of aver- 
age magnetisation in the dynamics 07j. From this fact it 
follows immediately that P+ = p. This result can be con- 
firmed by an explicit calculation. Starting with fixed mag- 
netisation xq = 2p— 1, the initial condition P(x, 0) = S(x— 
xq) broadens under the diffusive dynamics (|1 T|) and leaves 
always non-zero probability of ending in either of the pos- 
sible stationary states. We already noted that J x P(x, r) dx 



(XT 



0.01 



0.001 




Fig. 2. Probability of reaching the stationary state in time 
larger than r, for case II, q = 2, N = 2000. The values of 
initial fraction p of opinions +1 are 0.1 (+) 0.2 (x) and 0.7 
(©)• 



is independent of time under the dynamics 1)17(1. There- 
fore, the asymptotic state is the following combination of 
the eigenvectors (|20() with c = 

lim P(x,t) = ^-— ^Six + ^ + ^^Six-l) (33) 



and the probability of ending in the state of all +1 is 
therefore simply P + = p. 



4.3 Distribution of votes 

As already stressed in Sec. [3J equation i|17|) describes also 
the evolution of the distribution of votes in the case of q ^> 
1 parties. We will present an argument how our results 
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may explain the empirical data, suggesting the l/n law 
for the distribution of votes. 

As stressed in the discussion following Eq. (|23|l . the 
time-independent solutions of Eq. 1)17(1 can behave either 
as 1 + x or (1 + x)^ 1 in the limit x — > — 1. However, the 
latter case was excluded by the requirement of normal- 
isability of the probability density. On the other hand, 
relaxing the normalisability condition, the functions 



<t>02(x) = 



1 



l + x 
1 

l-x 



(34) 
(35) 



are solutions of (|22|) with eigenvalue c ~ 0. (Of course, any 
linear combination of them is also solution with c = 0). 

How should be any of these additional solutions inter- 
preted? The zero eigenvalue suggest that the function is 
stationary in time. However, it is not normalisable, there- 
fore this solution cannot be reached from any initial con- 
dition. But if the distribution P n (x,r) is close to 4>oi(%) 
(or 4>Q2{x)) in some interval / of x, it is probable that it 
P n (x,r) will remain close to (|34|) (or (|35|l . respectively) 
for certain period of time, while the interval / will grad- 
ually shrink and eventually disappear. Therefore, we may 
suggest <(3"4"1) and (|33|) as a metastable states, or long-lived 
transient states. 

This may explain the observation from simulations per- 
formed in 0] . In this work, the distribution of the type 1 /n 
is obtained in a suitably chosen transient regime, in cer- 
tain range of n. As x = In — 1, the behaviour of i|34|> at 
x — ► — 1 corresponds precisely to 1 jn behaviour for small 
n. 

A slightly more rigorous variant of the above argument 
is also possible. Imagine now that the political system rep- 
resented by the set of opinions 5" is not closed, but new 
opinions may appear, replacing other ones which vanish. 

Indeed, the current induced by the dynamics of case II 
can be read off from Eq. (|17|l 



7 



d_ 

dx 



[(1-x 2 )P(x,t)} 



(36) 



and by insertion of the solution l|34|) we deduce that there 
is homogeneous flow j = +1 outward the value x = — 1, 
i. e. n = 0. We may interpret this flow as a consequence 
of an external source placed somewhere close to the point 
x = — 1, i. e.n = 0. Such a source accounts for the influx 
of new opinions, or new parties, into the system. It is very 
reasonable to assume that the source is placed at very 
small values of n, as new subjects are likely to gain little 
support initially. 




Fig. 3. Average time of reaching the stationary state in dy- 
namics of case I, q = 2. The system size is N = 2000 (+) and 
N = 4000 (x). The line is the analytic prediction of Eq. (tT2t . 




Fig. 4. Relaxation time toward the stationary state in dy- 
namics of case I, q = 2. The system size is N = 2000 (+) and 
N — 4000 ( x ) . The horizontal line is the analytic prediction of 
Eq. 



focus was on the dynamical properties, namely the dis- 
tribution of times needed to reach the homogeneous sta- 
tionary state. We show in figures ^ and [21 the probabilities 
P^(t) that the time r s t to reach the stationary state is 
larger that r. We can clearly see that the probability de- 
cays exponentially with t in both cases I and II. 

Let us discuss the case I first. Following the analyti- 
cal expectation (|15|l we can fit the exponential tail of the 
distribution as 



5 Comparison with numerical simulations 

We performed numerical simulations of the Sznajd model 
on fully connected network according to algorithms de- 
scribed in sections 12.21 fcase I) and 12. 31 (case II). The main 



P>(r)~exp - 



'T'rclax 



(37) 



The results for (T s t) can be seen in Fig. EI compared with 
the analytical prediction of Eq. H12|l . Similarly in Fig.0|we 
can compare the fitted relaxation time with the analytical 
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Fig. 5. The fitted parameter to for reaching the stationary 
state in dynamics of case II, q — 2. The system size is N = 
2000. The line represents the formula 14011 . In the inset, the 
fitted first two relaxation times r r o and r r i are shown. The 
horizontal lines are corresponding analytical predictions from 
Eq. GU. 

result. Both (r s t) and r rc i ax agree satisfactorily with the 
analytical predictions. The deviations around the value 
p = 0.5 are due to finite size effects; the comparison of the 
results for system sizes N = 2000 and N = 4000 supports 
this interpretation. From Eq. I|12(l we can see that (T s t) 
diverges logarithmically for N — > oo. This is confirmed by 
the simulation data which fall onto single curve in Fig. [3] 
for different system sizes. 

Now let us turn to the case II. The equation (|29|) yields 
the leading term in the tail of the distribution (r) and 
in principle also the corrections to it. As the functions 
4> c (x) are odd for c = q with odd I, the next non-zero cor- 
rection will come from the eigenvalue C2 = 12. Therefore, 
we expect the behaviour 

P>(t) ~ exp ( ) + oi exp(- — ), r ^ oo (38) 

with 

T r0 = -, Trl = — . (39) 

As in the initial condition Po(^) = S(x — xo) we have 
xq = 2p — 1, we can deduce from Eq. H30[) the following 
estimate 

r ~ In y/6p(l-p) . (40) 

We can see from Fig. [S] that it corresponds well to the 
numerical data. In the inset of Fig. we can also see the 
fitted relaxation times r r o and t t %. Also here the corre- 
spondence with analytical prediction (|39(1 is good. 

6 Conclusions 

We formulated a mean-field version of the Sznajd model 
of opinion formation by putting it on a complete graph. 



Solving the underlying diffusion equations we found ana- 
lytical results for several dynamical properties, as well as 
exact long-time asymptotics. The results differ substan- 
tially in the two cases studied: first, the original Sznajd 
model, where a cluster of identical opinions is necessary to 
persuade others to join them, and second, the Ochrombel 
simplification, where also isolated agent can persuade oth- 
ers. Dynamical phase transition was found analytically in 
the original Sznajd model, while in the Ochrombel version 
it is absent. This finding agrees with previous numerical 
results. 

The approach to stationary state was the main con- 
cern of our calculations. We found that the distribution of 
times to reach the stationary state has an exponential tail 
which we were able to calculate analytically. In the case 
of Ochrombel simplification, wc obtained also the correc- 
tions and a formula which gives in principle the whole dis- 
tribution. We compared the analytical results for the tail 
(and in the Ochrombel case also for the first correction) 
with numerical simulations and we found good agreement. 
The method of adjoint equation enabled us to find ana- 
lytically the average time to reach the stationary state, in 
both cases. 

We found also another signature of the phase transition 
in the original Sznajd model, expressed by the divergence 
of the average time to reach the stationary state. Con- 
trary to the Ochrombel case, in the original Sznajd model 
the average time needed for reaching the stationary state 
blows up logarithmically with increasing system size. This 
finding was also confirmed in our numerical simulations. 

The analytical treatment provided an explanation of 
the 1/n distribution of votes, documented in Brazilian 
elections. We found that this behaviour corresponds to 
long-lived transient state of the dynamics of the Sznajd 
model with large number of possible opinions, or alterna- 
tively to the dynamics of an open version of the Sznajd 
model, where new opinions may continuously emerge. 

This work was supported by the project No. 202/01/1091 of 
the Grant Agency of the Czech Republic. 



Appendix A : Sznajd model on an arbitrary 
social network 

Our system is composed of N agents placed on nodes of a 
social network, represented by the graph A = (_T, E) where 
r is the set of nodes and E set of edges, i. e. unordered 
pairs of nodes. For a node i € F wc denote = {j £ 
r\(i,j) £ E} the set of neighbours of i. 

The opinion of the agent i s denoted a,-. The state of 
the system is described by the set of opinions of all the 
agents, S = [01, a2, <?n] G S r . The variable E(t) per- 
forms a discrete-time Markov process, defined as follows. 

In the case I we iterate the following three steps. First, 
choose i £ r at random. Then, choose j £ _Tj randomly 
among neighbours of i. If (7j(t) ^ ^ (i), nothing happens. 
However, if <7j(t) = <Jj(t), we will choose randomly one 
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of the common neighbours k G -Tj PI Tj \ and set 

(T k (t+1) =(Ji{t). 

In the case II we choose i <G r at random. Then, choose 
j G r"j randomly among neighbours and set crj(t + 1) = 
er<(t). 

If the graph is random and densely connected, we may 
approximate it by the complete graph with N nodes, i. e. 
for each pair of nodes i,j G P there is an edge connecting 
them, G E. It means that the set of neighbours of a 
node i E r is Ti = r \ {i}. This is a kind of a mean-field 
approximation. 



Appendix B : Finding the eigenvectors 

We can look for the solution of the equation (|22[l in the 
form of power series 

oo 
1=0 

and find the recurrence relation for the coefficients 



JZ+2 



1 



(i+m 



bi 



(B.2) 



We should distinguish two cases. Either the sequence 
of coefficients bi contains non-zero values for arbitrarily 
large I, or it is truncated at some order and (|B.1|) becomes 
a polynomial. In the former case the solution behaves as 
(j)c(x) ~ (1 — a; 2 ) -1 at x — » ±1 and must be excluded. The 
latter case is possible only if 

c = q = (Z + 1)(Z + 2) (B.3) 

for some I > 0. Moreover, in order to have a solution in 
the form of a polynomial, we require that b\ = if I in the 
equation 1|B.3|) is even, and bo = if I in the equation l|B.3(l 
is odd. The following table lists the solution for several 
lowest eigenvalues (taking bo = 1 for even I and b\ = 1 for 
odd I). 



Q 
2 
6 

12 
20 
30 



<f) c (x) 

1 

X 

1 - 5a; 2 
7 3 

1 - Ux 2 + 21x 4 



(B.4) 



In fact, up to a multiplicative constant, the functions 4> c {%) 
are Gegenbauer polynomials (531154) . 

The same procedure can be used for finding the eigen- 
vectors of the adjoint operator, solving Eq. I|25|l. We ex- 
pand the function tp c (x) in power series 



^c(z) =^2 d i x ' 



(B.5) 



and find the recurrence relation for the coefficients 

di+2 - Q+m+2) di ■ 



(B.6) 



Again we conclude that the only acceptable values of c are 
given by condition c = cj = (I + + 2) for I = 0,1, 2, ... 
and in this case the eigenvectors are polynomials of order 
1 + 2 in the variable x. The following table lists the solution 
for lowest eigenvalues (taking d a = 1 for even / and d\ = 1 
for odd I). 



1 


Q 




1p c (x) 





2 


1 - 


x 2 


1 


6 


x — 


x 3 


2 


12 


1 - 


6x 2 + 5x 4 


3 


20 


x — 


10 3 7 5 

t + r 


4 


30 


1 - 


15x 2 + 35a; 4 - 21a; 6 



(B.7) 



1=0 



It is important to note that the set of right eigenvalues 
coincides with the set of left eigenvalues, which proves 
consistency of our approach. 

Note that neither 4>c{%) nor ipc{z) are orthogonal poly- 
nomials. Instead, they are mutually orthogonal, i. e. 

(t> c (x)?p c i (x)dx = for c =/= d . This is due to the fact 
that the operator C is not self-adjoint. 
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